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^ ■ In the recent years we have seen a lot of activity around systems and 

Tjlj- I experiments, like DIS at HERA or the heavy-ion experiments at RHIC, 

i involving a large number of partons due to the high enery and/or the high 

CN| ' number of participants of those experiments. The main problem in this 

regime is that of the high parton densities. In fact, in most of the mod- 

' els of multiparticle production, two contributions to the multiplicity are 
considered: one proportional to the number of participant nucleons, Npart, 
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and a second one proportional to the number of collisions, N^l^^. In order 



to get the right multiplicities at RHIC it is necessary to lower the second 
contribution. A possible mechanism for this is the saturation. Here, I 
I am going to review the saturation of parton densities in the initial state, 

1^ ■ in two different frameworks: the Colour Glass Condensate and the string 

clustering. 



^ ; PACS numbers: 12.38.Mh, 24.85. -fp, 25.75.Nq 



1. THE COLOUR GLASS CONDENSATE 

1.1. Parton saturation at small x and the saturation momentum 

It has been much activity in the last years trying to understand the 
physics of nuclear and hadronic processes in the regime of very small Bjorken's 
X (very high energy). The main problem in this regime is that of the high 
parton densities. At high energy, the QCD cross sections are controlled by 
small longitudinal momentum gluons in the hadron wavefunction, whose 
density grows rapidly with increasing energy or decreasing x, due to the 
enhancement of radiative process. If one applies perturbation theory to this 
regime, one finds that, by resumming dominant radiative corrections at high 
energy, the BFKL equation leads to a gluon density that grows like a power 
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of s and in consequence to a cross section that violates the Froissart bound. 
Nevertheless, the use of perturbation theory to high-energy problems is not 
obvious. In fact, the BFKL and DGLAP equations are linear equations that 
neglet the interaction among the gluons. With increasing energy, recombi- 
nation effects -that are non-linear- favored by the high density of gluons 
should become more important and lead to an eventual saturation of parton 
densities. 

These effects become important when the interaction probability for the 
gluons becomes of order one. Taking "q^" as the transverse size of the 

gluon and ^'^^^^ ^ as the density of gluons, the interaction probability is 
expressed by 

OgNr xG(x,Q'^) 

Equivalently, for a given energy, saturation occurs for those gluons having a 
sufficiently large transverse size r'j_ larger than some critical value 

1/Qs{x, A). So the phenomenon of saturation introduces a characteristic 
momentum scale, the saturation momentum Qs{x,A), which is a measure 
of the density of the saturated gluons, and grows rapidly with 1/x and A (the 
atomic number). The probability of interaction -that can be understood as 
"overlapping" of the gluons in the transverse space- becomes of order one 
for those gluons with momenta ^ Qs{x, A) where 

2. xG{x,Ql) _ (colour charge)^ 

Qs{x,A) = asNc ^2 — = • (2) 

vr/t^ area 

For Q2 <Q'i^(x,A), the non-linear effects are essential, since they are ex- 
pected to soften the growth of the gluon distribution with r = ln(l/j;). For 
a nucleus, xGa{x, Q1) oc A and t\:R\ a A^/^, so eq. © predicts Ql oc A^/^. 
One can estimate the saturation scale by inserting the BFKL approximation 
into eq. . This gives (with (5 ~ 1/3 and A ~ ca.g in a first approximation) : 



Ql{x,A) ^ A^x-^, c= [ - /3 + ^/3(/? + 8u;) ]/2 = 4.84... , (3) 

which indicates that an efficient way to create a high-density environment 
is to combine large nucleai with moderately small values of x, as it is done 
at RHIC. 

This equation also shows that for sufficiently large energy or x small 
enough, Ql{x^ A) » ^}^CD ™^ Ois{Qs) !> which characterizes the regime 
of weak coupling QCD. But although the coupling is small, the effects of 
the interactions are amplified by the large gluon density: at saturation, 
Ga{x.,Q1) ~ l/as{Qs) S> 1, so the gluon modes have large occupation 
numbers, of order (corresponding to strong classical fields A ~ 

which suggests the use of semi-classical methods. 
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1.2. The effective theory for the CGC: the Renormalization Group 

Equation. 

One can write a classical effective theory based on this general idea: the 
fast partons -valence quarks with large longitudinal momentum- are consid- 
ered as a classical source that emits soft gluons -with smaller longitudinal 
momenta- which are treated as classical colour fields A[p]. The fast partons 
move nearly at the speed of light in the positive x"*" direction, and generate 
a colour current J'^ = 5^+/?. By Lorentz contraction, the support of the 
charge density p is concentrated near the light-cone longitudinal coordinate 
x~ = 0. By Lorentz time dilation, p is independent of the light-cone time 

The Yang Mills equation describing the soft gluon dynamics reads 

j^^puf. ^ 5/^+p(^-,x) . (4) 

Physical quantities, as the unintegrated gluon distribution, are obtained as 
an average over p: 

< A^{X) A^{Y) >, = J D[p] W,[p] A^[p]{X) A'[p]{Y) , (5) 

where yl*(X) corresponds to the classical solution for a given p, and 
is a gauge-invariant weight function for p. What we are doing is a kind 
of Born-Oppenhaimer approximation: first, we study the dynamics of the 
classical fields (Weizsacher- William fields) for a given configuration p of the 
colour charges, and second, we average over all possible configurations. 
For the classical solution we find: 

F+\x-,x^) = 5{x~)^Vix^){d'V{x^)^) = d+A' , (6) 

where V{x_i_) is the Wilson line 

V\x±) = Peyip\^igJdx~A+{x-,x)} (7) 

and ^^[/o] is the solution of the equation of motion Q) in the covariante 
gauge: — V^^^"*" = p. 

The weight function VFr[p] is obtained by integrating out the fast par- 
tons, so it depends upon the rapidity scale r = ln(l/3;) at which one con- 
siders the effective theory. This can be taken into account via a one-loop 
background field calculation, and leads to a renormalization group equa- 
tion (REG) for VFt-[/9] which shows how the correlations of p change with 
increasing r. Schematically: 

dWr[p] 



dr 
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This is a functional diffusion equation, where the kernel x[p\ plays the role 
of the diffusion coefficient in the functional space spanned by p{x~,x±). 
This kernel is positive definite and non-linear in p to all orders. It depends 
upon p via the Wilson line (jT)). 

As it has been said above, the physical quantities, as the gluon density, 
are obtained as an average over p: 

n(x, k^) ^ ^ oc (F+'(fc^)F+'(-fe^)). . (9) 

In order to calculate this average, we need to solve the REG ((HI). Approxi- 
mate solutions to this equation can be obtained in two limiting cases: 

• At low energy, or large transverse momenta ^ Q'^{x), we are in 
a dilute regime where fields and sources are weak, and the Wilson lines 
can be expanded to lowest order in V^x_i) !^ I + ig J dx'^ A~^{x~ ,x±). 
In this case, the REG equation reduces to the BFKL equation, and the 
gluon density of eq. @ grows both with and 1/x (Bremsstrahlung): 
n{x,k±) ~ 

• At high energies, or low momenta k'^ ^ Q'^{x), the colour fields are 
strong, ^4+ ~ 1/g, so the Wilson lines rapidly oscillate and average away to 
zero: F ^ yt ~ 0. Then the kernel X becomes independent of p, and we 
obtained a gluon density that increases linearly with the evolution "time" 
T = ln{l/x): n{x,k±) ~ ^ In ^^^^ oc In^. That is, we find saiuraiion for 

the gluon density, that grows logaritmically with the energy since r ~ Ins: 
unitarity is restored. 

We call the high density gluonic matter at small-x described by this 
effective theory a Colour Glass Condensate (CGC) P: Colour since gluons 
carry colour under SU{Nc); Glass since we have a random distribution of 
time-independent colour charges which is averaged over in the calculation 
of physical observables, in order to have a gauge independent formulation; 
and Condensate because at saturation the gluon density is of order l/og, 
typical of condensates, so we have a system of saturate gluons that is a Bose 
condensate. 



1.3. Phenomenology at RHIC 

The particle production in RHIC collisions has been analyzed from the 
perspective of the CGC, considering it as a pertinent description of the 
initial conditions. Taking into account that the multiplicity is proportional 
to the number of gluons -parton-hadron duality-, the centrality dependence 
of multiparticle production has been related to the density of gluons |2], 
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that at saturation (see eq. ((^): 

where t^R\ oc corresponds to the nuclear overlap area, and Qg is the 

1 /3 

saturation momentum for the considered centrality, Ql{x,A) oc Np^^.^. To 
compute the centrality dependence, it is necessary to know the evolution 
of the gluon structure function, wich is governed by the DGLAP equation. 
Taking l/as(Qs) ~ ^^^(QsMqcd) ~ ^T^^part, one finally finds that the 
multiplicity per participant behaves as In A'pa^^. 

Besides, it has been obtained from saturation a proportionality between 
the mean transverse momentum and the multiplicities, 

2 1 diV 
<^^>'^^d^- ^''^ 

This observation indicates that the px broadening seen in elementary and 
heavy ion collisions results from the same physics, the intrinsic generated 
PT broadening in the partonic phase [S]. 



2. STRING CLUSTERING 

In the clustering approach, the colour strings created in the nuclear 
collisions are considered as effective sources with a fixed transverse area, 
r_|_ ~ 0.2 fm. Notice that this radius coincides with the estimate satura- 
tion momentum of the CGC at RHIC. If the string density is high enough, 
the strings overlap, forming clusters very much like disks in continuum 
two-dimensional percolation theory. In order to calculate the physical ob- 
servables, as the particle multiplicity or the mean transverse momentum, 
we need to study the dynamics of those clusters. 

We assume that a cluster of n strings behaves as a single string with 
a higher colour field Qn, corresponding to the vectorial sum of the colour 
charges of each individual Qi string. The resulting colour field covers the 
area Sn of the cluster. As = Qi)'^i and the individual string colours 
may be oriented in an arbitrary manner, the average QuQij is zero, so 
Qn — '^Qi- Qn depends also on the area Si of each individual string that 
comes into the cluster, as well as on the total area of the cluster Sn, Qn = 
\J^^§^Qi- We take Si constant and equal to a disc of radius r±_ ~ 0.2 fm. 
Sn corresponds to the total area occupied by n discs^. Knowing the colour 

^ Notice that if the strings don't overlap, Sn = nS\ and Qn = nQ\, so the strings act 
independently. On the contrary, if they fully overlap, Sn = Si and Qn = \fnQ\. 
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charge Qn, one can compute the multipHcity ^„ and the mean transverse 
momentum < p\ of the particles produced by a cluster of n strings. One 
finds 

Atn = W^/^l, < Pt >n= d — < Pt >i (12) 

where fii and < p^ >i ^-re the mean multiplicity and mean pj, of particles 
produced by a simple string. In the saturation limit, i. e. all the strings 
overlap into a single cluster that occupies the whole interaction area, one 
gets the following scaling law that relates the mean transverse momentum 
and the multiplicity per unit rapidity and unit transverse area: 

/ 2 ^ Si <Pt>i .,on 

< Pt >aa= t; I^AA • (13) 

Saa fJ-i 

This scaling relation is similar to the one obtained in the framework of the 
CGC when the initial gluon density saturates. 

Moreover, in the limit of high density 77 = NgSi/ Saa, one obtains 

nSi rj _ 1 

^ 57^" l-exp(-7?) = 
and the equations ()12() transform into the analytic ones [H] 

/i = NsF{r])fii , < pf. >= < >i (15) 

where fi and < Px > are the total multiplicity and mean momentum, and 
Ng is the total number of strings created in the collision. 

In order to study the transverse momentum distribution, one needs the 
distribution f{x,mT) for each string or cluster, and the cluster size distri- 
bution W{x). For f{x,mT) we assume the Schwinger formula, /(x,mr) = 
exp(— mj.xF(?7)), used also for the fragmentation of a Lund string. In this 
formula x is related to the string tension, or equivalently to the mean trans- 
verse size of the string. The weight function W{x) obeys a gamma distri- 
bution 6 . Assuming that a cluster behaves similarly to a single string but 
with different string tension, that depends on the number of strings that 
come into the cluster, we can write for the total pT distribution 

/{mr) = J W{x) f{x,mT), W{x) = ^(^x)^"! exp(-7x) . (16) 

Performing the integral in eq. ()16|) we obtain 

diV _ dA^A:-l 1 
d^ " d^^^'^^i I ' ^ ^ 

where k can be obtained from the fluctuations in multiplicy. 
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3. SOME SIMILARITIES: The results from string clustering 

and from the CGC 

To finish, let us remember some of the similarities of the explained ap- 
proaches: 

• In the clustering approach, when taking the saturation limit -all the 
strings overlap into a single cluster that occupies the whole nuclear overlap 
area-, one finds that the particle multiplicity of a central collision, haa, 

behaves as fj,AA = fJ-n = \J "^^/^i = \/^^r^^i- Taking into account that 
the number of strings produced in the nuclear collision, Ng^ is proportional 
to the number of inelastic nucleon-nucleon collisions, Nf-oU ~ A^^^, and 
Saa corresponds to the nuclear overlap area, Saa ~ A^^^, we obtain a 
multiplicity that scales with the number of participants A. This coincides 
with the multiplicity obtained in first approximation - without evolution - 
in the framework of the CGC (eq. (fTfl|) ). 

• In the CGC, there is a relation between the mean transverse mo- 
mentum and the multiplicity, which is developed in eq. (jllj) . This relation 
shows that at saturation the mean transverse momentum should scale with 
the multiplicity per unit rapidity and unit transverse area. This coincides 
with the proportionality relation obtained in the clustering model (eq. H13() ). 

• In both approaches, the initial state interactions -gluon saturation in 
the CGC or clustering of strings- produce a suppression of high pT and 
multiplicities. On the contrary, in the framework of the jet quenching phe- 
nomena, the energy loss of the jet with a hot and dense medium produces 
additional soft gluons that would fragment into hadrons increasing the mul- 
tiplicities, unless strong shadowing occurs in the gluon structure functions. 
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